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The order parameter of a condensate with two internal 
states can continuously distort in such a way as to remove 
twists that have been imposed along its length. We observe 
this effect experimentally in the collapse and recurrence of 
Rabi oscillations in a magnetically trapped, two-component 
Bose-Einstein condensate of 8r Rb. 



Quantization and persistence of current in supercon- 
ductors and superfluids can be understood in terms of 
the topology of the order parameter. Current arises from 
a gradient in the phase of the order parameter. Quantiza- 
tion of flow around a closed path is a consequence of the 
requirement that the order parameter be single- valued; 
metastability, or "persistence", arises from the fact the 
number of phase windings (the multiple of 2n by which 
the phase changes) around the path can be changed only 
by forcing the amplitude of the order parameter to zero 
at some point. If the energy this requires exceeds that 
available from thermal excitations, then the current will 
be immune to viscous damping. This familiar argument 
relies, however, on the order parameter's belonging to a 
very simple rotation group. The order parameter in su- 
perfluid 4 He, for example, is a single complex number. 
Its phase can be thought of as a point lying somewhere 
on a circle, and is subject to the topological constraints 
mentioned above. The order parameter of a more com- 
plicated superfluid, on the other hand, will in general be 
capable of ridding itself of unwanted kinetic energy by 
moving continuously through a higher-dimensional order- 
parameter space in such a way as to reduce, even to zero, 
its winding number. Presumably this ability will reduce 
a superfluid's critical velocity; in the limit that the order- 
parameter space is fully symmetric, the critical velocity 
may even vanish [jjj. 

In this paper, we discuss experiments on a gas-phase 
Bose-Einstein condensate with two internal levels 0. 
This is equivalent to a spin- 1/2 fluid: the order param- 
eter has SU(2) rotation properties. A differential torque 
across the sample is applied to the order parameter so 
that with time it becomes increasingly twisted. Eventu- 
ally the sample distorts through SU(2) space so that the 
steadily applied torque now has the effect of untwisting 
the order parameter, which returns nearly to its unper- 
turbed condition. The pattern of twisting and then un- 



twisting is manifested experimentally as a washing-out 
followed by a recurrence of an extended series of oscilla- 
tions in the population between the spin states. Related 
behavior has been previously observed in A-phase 3 He 
p|; a major difference in this work is that we can di- 
rectly observe components of the order parameter with 
temporal and spatial resolution. 

Magnetically confined 87 Rb can exist in a superposi- 
tion of two internal states, known as |1) and |2) g. The 
two internal states are separated by the relatively large 

Rb hyperfine energy, but in the presence of a near- 
resonant coupling field the states appear, in the rotating 
frame, to be nearly degenerate. The condensate can then 
dynamically convert between internal states. The order 
parameter for the condensate is the pair of complex field 
amplitudes <&i and $2 of states |1) and |2). Evolution 
of these fields is governed by a pair of coupled Gross- 
Pitaevskii equations which model the coupling drive, the 
external confining potential, kinetic energy effects and 
mean-field interactions p~n|. The SU(2) nature of the 
order parameter (<&i,$2) is more evident if we write 



$1 = cos(6»/2)e 



and 



$2 = sin(0/2)e 
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where 0, 0, n t , and a are purely real functions of space 
and time. and <j> give the relative amplitude and phase 
of the two internal components, and may be thought of 
respectively as the polar and azimuthal angles of a vector 
whose tip lies on a sphere in SU(2) space. The total 
density and mean phase, n t and a respectively, remain 
relatively constant p| during the condensate evolution 
described in this paper. 

The apparatus has been previously described |5|,||,[|. 
The starting point for the measurements is a magneti- 
cally confined cloud of ~ 8 x 10 5 evaporatively cooled, 
Bose-Einstein-condensed 87 Rb atoms near zero temper- 
ature. The combined gravitational and magnetic poten- 
tials pQ l yield an axially symmetric, harmonic confining 
potential V\ (V2) for particles in the |1) (|2)) state, in 
which the aspect ratio of the axial oscillation frequency 
in the trap to the radial frequency u) z /u) r can be varied 
from 2.8 to 0.95 pfl. V\ and V% are nearly identical but 
can optionally be spatially offset a distance z$ in the axial 
direction p2J. The coupling field has a detuning S from 



the local |1) to 2} resonance. If zq is nonzero, 6 depends 
on the axial position z, with S(z) — S(z = 0) linear in z 
and in zq. The strength, characterized by the Rabi fre- 
quency f2, of the coupling field also varies with an axial 
gradient fll3| . 

We are able to measure the population of both spin 
states nondestructively using phase-contrast microscopy 
jnLfLq]. We tune the probe laser between the resonant 
optical frequencies for the |1) and |2) states. Since the 
probe detuning has opposite sign for the two states, the 
resulting phase shift imposed on the probe light has op- 
posite sign, such that the |1) atoms appear white and 
the 1 2) atoms appear black against a gray background on 
the CCD array. We can acquire multiple, nondestruc- 
tive images of the spatial distribution of the |1) and |2) 
atoms at various discrete moments in time, or we can 
acquire a quasi-continuous time record (streak image) of 
the difference of the populations in the |1) and |2) states, 
integrated across the spatial extent of the cloud. 

The effect of the coupling drive is to induce a preces- 
sion of the order parameter at the local effective Rabi 
frequency Q c s(z) = (£!(z) 2 + 8{z) 2 ) 1 / 2 . In a prelimi- 
nary experiment, we chose parameters so as to make $7 c ff 
nearly uniform, with lo z — 2n x 63Hz, to r — 2tt x 23 Hz, 
il ~ 2tt x 340 Hz and 6(z) ~ 0. A condensate at near-zero 
temperature was prepared in the pure |1) state. The cou- 
pling drive was then turned on suddenly, inducing an ex- 
tended series of oscillations of the total population from 
the |1) to the |2) state ("Rabi oscillations") [Fig. §. The 
robustness of the Rabi oscillations is proof that our imag- 
ing does not significantly perturb the quantum phase of 
the sample |T(| (population transfer via Rabi oscillations 
is phase-sensitive). 

If there is an axial gradient to fl e g, then a relative 
torque is applied to the order parameter across the con- 
densate, which can cause a twist to develop along the 
axial direction. If we naively model the sample as a col- 
lection of individual atoms, each held fixed at its respec- 
tive location, then the order parameter at each point in 
space rotates independently at the local effective Rabi 
frequency, f2 e ff- In Fig- |2|( a ) we see the implications 
of the "fixed-atom" model for Q = 2n x 700 Hz and 
S(z) — 2-7T x (100 + 14z)Hz with z in microns: a twist 
develops in the order parameter which leads to a wash- 
ing out of the Rabi oscillations [Fig. 0(c)]. In contrast, 
the kinetic energy provides stiffness for a true conden- 
sate. Simulation of the condensate |L7| shows that, for 
early times, the order parameter begins twisting as in the 
"fixed-atom" model but the twisting process self-limits 
about 40 ms into the simulation. At this point there is 
nearly a full winding across the condensate. Thereafter, 
though the two ends of the order parameter continue to 
twist with respect to one another, the order parameter 
has been sufficiently wrapped around the SU(2) sphere 
that the effect of further torque is to return the conden- 
sate close to its unperturbed condition [Fig. ||(d)]. The 



Rabi oscillations exhibit a corresponding revival [Fig. 
0(e)]. The factor driving the untwisting process is the in- 
creasing kinetic energy cost associated with an increasing 
twist in the order parameter. For a simple U(l) order pa- 
rameter, continuously increasing the winding ultimately 
results in a "snap" , in which the order parameter is driven 
to zero and a discontinuous (and presumably dissipative 
fll8|) process releases the excess windings. The revival in 
the present case is made possible by the larger rotation 
space available to a two-component cloud. 

Under experimental conditions similar to those of the 
simulations in Fig. 0, we have observed as many as 
three complete cycles of Rabi-oscillation decay and re- 
currence. These data appear in Ref. |p_9[] . In this paper 
we present data that correspond to the case of a more 
vigorous twisting. We increase the axial dimension of 
the condensate cloud by a factor of four; the kinetic en- 
ergy cost of twisting the condensate is correspondingly 
lower, so that at the point in time when the condensate 
is maximally distorted there are four windings across the 
cloud. The parameters of the experiment were as follows: 
U) r = <jJ z — 2tt x 7.8 Hz and mean tt c s = 2ir x 225Hz. 
There was a gradient in both S and f2 across the 54/im 
axial extent of the cloud, resulting in a ~ 2ir x 60Hz 
difference in n e g from top to bottom of the condensate. 
The result of the experiment is seen in Fig. 0. The ob- 
served recurrence of the Rabi oscillations at 180 ms [Fig. 
0(a)], when corrected for overall decay of the cloud, cor- 
responds to 60 percent contrast. We find it remarkable 
that the distorted order-parameter field seen in Fig. 0(b) 
at times 65 and 75 ms should find its own way back to a 
nearly uniform configuration. 

The simulations qualitatively reproduce the integrated 
number and state-specific density distributions observed 
in the experiments. For large inhomogeneity in Q c g , how- 
ever, the simulations predict the development of small- 
scale spatial structure not observed in the experiment. 
The simulations contain no dissipation, whereas finite- 
temperature damping may occur in the experiment |20| . 

Heuristically, what value do we expect for the recur- 
rence time irccur for the data in Fig. g? The differ- 
ence in Q e ff from the top to the bottom of the con- 
densate is about 60 Hz. From the data in Fig. ||(b) 
we see that the recurrence occurs only after four wind- 
ings have one-by-one been twisted in and then twisted 
out of the condensate. A rough estimate then would be 
irccur = (4 + 4)/60Hz= 133 ms, shorter than the observed 
value of 180 ms, but reasonable given that edge effects 
have been neglected. 

An interesting theoretical challenge would be to de- 
velop simple arguments that would allow an a priori 
prediction of the spacing of the windings at the instant 
of maximum twist. For particularly strong torques, one 
might expect the total density to be suppressed to zero 
along a plane transverse through the cloud, so as to al- 
low for discontinuous relaxation of the order parameter. 
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Indeed we have seen such behavior in numerical Simula- [6] 

tions. Under what conditions should suppression of total 
density, rather than continuous evolution through SU(2) I 7 ] 

space, be the preferred mode of relieving accumulated 
stress? 

We have observed that the presence of the coupling 
drive need not result in population oscillations. For any 
given frequency and coupling strength, there are two 
steady-state solutions which are completely analogous 
to the dressed-states solutions of the single-atom prob- 
lem pi]]. We have been able experimentally to put the 
condensate in such states via an adiabatic process: the 
strength of the drive is increased gradually from zero, and 
the frequency (initially far detuned) is gradually ramped 
onto resonance. The resulting "dressed condensate" is 
extremely stable — after the ramp is complete, the cloud 
remains motionless in a near-equal superposition of the 
two bare states. The twisting experiment discussed in 
this paper can be thought of as the evolution of a highly 
nonlinear superposition of the two dressed states. It 
would be interesting to work in the opposite limit and 
explore the spectrum of sma/Z-amplitude excitations on a 
dressed state. Analogies to (i) the transverse zero-sound 
in He 3 B-phase p2| a nd (ii) spin-waves in spin-polarized 
atomic hydrogen p3| should also be studied [ p4[ . 

The ability to fundamentally alter the topological 
properties of a condensate has already proven useful. 
When the coupling drive is turned off, the SU(2) proper- 
ties vanish, and states |1) and |2) become distinct species, 
each forced to live in its separate U(l) space. Figure || 
illustrates how this can be used to change $i(z) in a con- 
trolled manner. With a variation p5j on this technique 
we have created a vortex-state condensate [E6[. 

We acknowledge funding from the ONR and the NSF, 
and useful conversations with Jason Ho and Seamus 
Davis. 



[14] 

[15. 
[16. 



[17. 

[is; 

[19 

[20 

[21 

[22 

* Quantum Physics Division, National Institute of Stan- [23 
dards and Technology. 

[1] P. Bhattacharyya, T. Ho, and N. Mermin, Phys. Rev. [24 
Lett. 39, 1290 (1977);T.-L. Ho, Phys. Rev. Lett. 48, 1837 
(1982). 

[2] For other recent experimental work on multi-component [25 
BEC, see J. Stenger et al, Nature 396, 345, (1998), and [26 
H.-J. Miesner et al, Phys. Rev. Lett. 82, 2228 (1999). 

[3] D. Paulson, M. Krusius, and J. Wheatley, Phys. Rev. 
Lett. 37, 599 (1976). 

[4] The states |1) and |2) correspond to the two zeeman sub- 
levels, F — l,m — — 1, and F — 2, m — 1 respectively. 
For further details on the system, see Ref. M. 

[5] E. A. Cornell, D. S. Hall, M. R. Matthews, and C. E. 
Wieman, J. Low Temp. Phys. 113, 151 (1998). 



D. S. Hall, M. R. Matthews, C. E. Wieman, and E. A. 

Cornell, Phys. Rev. Lett. 81, 1543 (1998). 

R. Dum, J. I. Cirac, M. Lewenstein, and P. Zoller, Phys. 

Rev. Lett. 80, 2972 (1998); R. J. Ballagh, K. Burnett, 

and T. F. Scott, Phys. Rev. Lett. 78, 1607 (1997); K.-P. 

Marzlin, W. Zhang, and E. M. Wright, Phys. Rev. Lett. 

79, 4728 (1997); J. Williams et al, Phys. Rev. A 57, 

2030 (1998). 

D. S. Hall et al, Phys. Rev. Lett. 81, 1539 (1998). 

M. R. Matthews et al, Phys. Rev. Lett. 81, 243 (1998). 

W. Petrich, M. H. Anderson, J. R. Ensher, and E. A. 

Cornell, Phys. Rev. Lett. 74, 3352 (1995). 

JR. Ensher Ph D, thesis, University of Colorado (1998). 

D.S. Hall et al." 



ond-mat/990345E 



The microwave coupling drive is actually accomplished 
by a two-photon transition. The quoted values of il and 
S refer to the combined effect of a two-frequency drive. 
See Ref. M. Since fi depends on the intermediate-state 
detuning, it can in general have a spatial dependence due 
to Zeeman shifts. £1 also depends on the (poorly charac- 
terized) polarization of the coupling fields, whose projec- 
tion onto the local quantization axis varies weakly across 
the cloud. In our comparisons with theory, the gradient 
in Q is an adjustable parameter. 

F. Zernike, in Nobel Lectures, Physics (Elsevier, Amster- 
dam, 1964), p. 239, text of 1953 Nobel Prize lecture. 
M. R. Andrews et al, Science 273, 84 (1996). 
The probe beam does impose a small, reproducible shift 
on the energy difference between the |1) and |2) states, a 
measurable shift for which we correct. 
The simulation is a numerical integration of the coupled 
Gross-Pitaevski equations in ||. See |19| for more detail. 
K. Schwab, N. Bruckner, and R. Packard, 386, 585 
(1997); O. Avenel, P. Hakonen, and E. Varoqaux, Phys. 
Rev. Lett. 78, 3602 (1997). 

J. Willia ms, R. Walser, J. C ooper, E.A. Cornell and M. 
cond-mat/9904399[ 

G. V. Shlyapnikov, and J. T. M. Wal- 
raven, Phys. Rev. Lett. 80, 2269 (1998); L. P. Pitaevskii 
and S. Stringari, Phys. Lett. A 235, 398 (1997); W. V. 
Liu, Phys. Rev. Lett. 79, 4056 (1997). 
P. B. Blakie, R. J. Ballagh and C. W. Gardiner, cond- 



Holland 

P. O. Fedichcv 



mat/9902110 



Y. Lee, T.M. Haard, J.B. Kycia, J. A. Sauls and W.P. 

Halperin, Nature 400, 431, (1999). 

N. Bigelow, J. Freed, and D. Lee, Phys. Rev. Lett. 63, 

1609 (1989). 

Obviously, oscillatory behaviour in coupled spin systems 

is a recurring theme. At present it is unclear to us how 

deep the analogy runs. 

J.E. Williams and M.J. Holland, private comm. 

MR. Matthews, B.P. Anderson, P.C. Haljan, D.S . Hall, 



C.E. Wieman, and E.A. Cornell, :ond-mat/9908209 




50 100 

Time (ms) 



FIG. 1. (a) With the trap parameters adjusted for high 
spatial uniformity in f2 e ff , we drive the coupling transition and 
record a streak-camera image of 60 Rabi oscillations between 
the |1) (white) and |2) (black) states. The vertical dimension 
of the figure is 80 /im. (b) The value of SN, the total number 
of atoms in |2) minus the total in |1), is extracted from the 
image in part (a). The contrast ratio remains near unity 
— observed loss of signal is due to overall shrinkage of the 
condensate through collisional decay. 
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FIG. 2. (a) We represent the polar-vector order parameter 
as an arrow in these simulations. The angle 6 from the vertical 
axis determines the relative population, and the azimuthal 
angle cj> is the relative phase of states |1) and |2) (see Eqs. 1). 
Each column in the arrow-array is at fixed time, and each row 
at fixed axial location. Cl is perpendicular to the plane of the 
page, so that a uniform, on-resonance Rabi oscillation would 
correspond to all the arrows rotating in unison, in the plane 
of the image. The tips of all the arrows are (on the relatively 
fast time scale of fi e ff) tracing out circles nearly parallel to 
the plane of the page (in our rotating-frame representation, 
small excursions out of the page are a consequence of finite 
detuning). In the figures, we "strobe" the motion just as 
the central arrow approaches vertical, to emphasize the more 
slowly evolving "textural" behavior, (b) The total density 
of the condensate nt maintains a Thomas- Fermi distribution 
(integrated through one dimension, as imaged) and changes 
only slightly during the evolution of the cloud, (c) In a simple 
model of individual, fixed atoms, a continuous inhomogeneity 
in f2 e ff will cause the Rabi oscillations in SN to wash out. 
(d) When a condensate is simulated [19], the kinetic energy 
causes the order parameter to precess through the full SU(2) 
space, coming out of the page to cast off the winding and thus 
reduce its kinetic energy, (e) The corresponding plot of SN 
shows that when the arrows are once more aligned, the Rabi 
oscillations recur. 
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FIG. 3. A condensate with large axial extent undergoes 
twisting, (a) The streak camera data shows a rapid decay in 
the Rabi oscillation in the integrated population difference, 
from full contrast at t=0 to near zero by t — 20ms. The 
oscillations recur at 180 ms. (b) Individual phase-contrast 
images (at distinct moments in time) of the spatial distribu- 
tion of |l)-state atoms show that the spatially inhomogeneous 
Rabi frequency is twisting the order parameter, cranking suc- 
cessively more windings into the condensate until by ~ 75 
ms four distinct windings are visible. Further evolution re- 
sults not in more but in fewer windings until, at time 180 ms, 
the order parameter is once more uniform across the cloud. 
Each image block is 100 /tm on a side, and the probe laser is 
tuned much closer to the |1) state than to the |2) state, (c) 
The numerical simulation reproduces the qualitative features 
of the corresponding experimental plot (a). The simulation 
used S(z) — 0, Q(z = 0) = 27r x 225Hz and a 2n x 60Hz spread 
in Q, across the extent of the condensate. 



